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Abstract 

We study Landau-Zener transitions in a dissipative environment by means of the quasiadiabatic propa- 
gator path-integral scheme. It allows to obtain numerically exact results for the full range of the involved 
parameters. We discover a nonmonotonic dependence of the Landau-Zener transition probability on the 
sweep velocity which is explained in terms of a simple physical picture. This feature results from a nontriv- 
ial competition between relaxation processes and the external sweep and is not captured by perturbative 
approaches. In addition to the Landau-Zener transition probability, we study the excitation survival prob- 
ability and also provide a qualitative understanding of the involved competition of time scales. 
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1. Introduction 

The transition dynamics in a quantum two-level 
system with a time-dependent Hamiltonian vary- 
ing such that the energy separation of the two di- 
abatic states is a linear function of time is a cen- 
tral problem since the early days of quantum me- 
chanics. It is commonly denoted as the Landau- 
Zener (LZ) problem, although it has been solved 
independently by Landau [l[ , Zener Q , Stiickelberg 
Q and Majorana Q in 1932 (for a more detailed 
discussion of the differences between the four ap- 
proaches, we refer the interested reader to a paper 
by Di Giacomo and Nikitin Q). Nonadiabatic tran- 
sitions at avoided level crossings are at the heart 
of many dynamical processes throughout physics 
and chemistry. They have been extensively stud- 
ied both theoretically and experimentally in, e.g., 
spin flips in nanomagnets @, solid state artificial 
atoms 0,0,0! nanocircuit QED [nj EH, adiabatic 
quantum computation fl2| , the dynamics of chemi- 
cal reactions [13j , and in Bose-Einstein condensates 
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in optical lattices |14j . Also the noncquilibrium dy- 
namics of glasses at low temperatures is dominated 
by the transition dynamics of avoided level cross- 
ings HE El, El Ell- 



in the pure Landau-Zener problem, two quantum 
states interact by a constant tunneling matrix ele- 
ment Ao. A control parameter is swept through 
the avoided level crossing at a constant velocity 
v, such that the energy gap between the two dia- 
batic states depends linearly on time. The Landau- 
Zener problem addresses the case when the system 
starts in the lower energy eigenstate in the infi- 
nite past and asks for the probability of finding the 
system in the lower energy eigenstate in the infi- 
nite future (a Landau-Zener transition). Certainly, 
for infinitely slow variation of the energy difference 
v — > 0, the adiabatic theorem states that no tran- 
sition between energy eigenstates will occur, since 
at any moment of time, the system will always be 
in an instantaneous eigenstate of the Hamiltonian. 
For v ^ 0, the probability Pq for no transitio n Il9 1 
is described by the Landau-Zener formula [H, El Elli I 
P (v,A ) =l-exp[-7rAg/(2t;)]. 

Most experimental investigations, on the one 
hand, actually deal with more complex systems 
than two-state systems and, on the other hand, do 



Preprint submitted to Elsevier 



March 3, 2013 



not necessarily start in the ground state. In close 
proximity around the avoided crossing in most cases 
only the two crossing states are important to de- 
scribe the dynamics. However, in general in such 
an approximate system we also have to deal with 
the system being initially in the excited state. In 
the pure Landau-Zener problem the probability to 
end up in the ground (excited) state when starting 
in the ground (excited) state are identical due to 
symmetry. 

However, in any physical realization, a quantum 
system is influenced by its environment leading to 
relaxation and phase decoherence during time evo- 
lution 2(J 21 1 . At low temperatures when the en- 
vironmental fluctuations are not thermally occu- 
pied, the Landau-Zener probability P (to end up 
in the ground state when starting in the ground 
state) should hardly be influenced since no phonons 
are available for the system to be in the excited 
state at asymptotically long times. However, spon- 
taneous emission is possible even at lowest temper- 
atures and the excitation survival probability Q (to 
end up in the excited state when starting in the 
excited state) is expected to decrease when the sys- 
tem has the time to relax during the driving. Thus 
especially at low driving speeds, the excitation sur- 
vival probability will be reduced since the system 
will decay, while the pure Landau-Zener mechanism 
predicts full survival of excitations when driving 
through an avoided crossing at low speeds. 

The dissipative Landau-Zener problem received 
a lot of attention 0, [H, [U [H, ^ [H, [H, Eil in 
the past 25 years due to its relevance for controlled 
quantum state preparation which became experi- 
mentally feasible in many physical realizations. Al- 
though the full problem is analytically unsolved, 
many limiting cases are analytically tractable and 
numerical approaches are available. Usually a dis- 
sipative environment causes fluctuations of the en- 
ergies of the diabatic states (longitudinal or diago- 
nal system-bath coupling) but occasionally the en- 
vironments can also cause transitions between the 
diabatic states (transversal system bath coupling). 
In this work, we focus on the more common case of 
longitudinal coupling. The transversal coupling has 
been treated in Refs. 24, 2f| for zero temperature. 



In the limit of very fast sweeps (nonadiabatic 
driving), a thermal heat bath has been shown not 
to influence the Landau-Zener probability. This fol- 
lows consistently from all approaches (for the cor- 
respondin g lim iting conditions of the other free pa- 
rameters) in m im m, m m m, sn . 



At low temperatures, the Landau-Zener (or, 
cquivalcntly, ground state survival) probability P 
(to end up in the ground state when starting in the 
ground state) is only influenced little by a diago- 
nally coupled bath. This was first shown by Ao and 
Rammer j24j. Kayanuma 22, 23 1 showed before 
that in the limit of slow fluctuations of the diabatic 
energies (as might be caused by a diagonally cou- 
pled bath at low temperatures), the Landau-Zener 
probability is not modified. At zero temperature, a 
diagonally coupled bath has strictly no influence on 
the Landau-Zener probability, as has been shown by 
Wubs et al. [25|. This is true not only for bosonic 
but also for spin environments (26]. In contrast, 
transversely coupled baths lower the Landau-Zener 
probability even at zero temperatures [25j and their 
effects depend on the details of the bath character- 
istics m. 

The strong damping limit can be reached in 
two ways. For large coupling between system and 
bath and finite temperatures, Ao and Rammer 24 1 
have shown that again a diagonally coupled bath 
does not influence the Landau-Zener probability. 
Conversely, in the high temperature limit under 
adiabatic conditions (i.e., slow sweeps), the two 
states are driven to equal population, P = PgD = 
|(1 — exp(7rAn/t>)). This has first been shown by 
Kayanuma [2z 



23| for large and fast fluctuations 
and later by Ao and Rammer (24j. Pokrovsky and 
Sun [28| have extended Kayanuma's result to trans- 
versely coupled environments. 

For the experimentally important parameter 
range of slow (adiabatic) sweeps and intermedi- 
ate temperatures but only weakly coupled envi- 
ronments, Pokrovsky and Sun (28l |. Kayanuma and 
Nakayama 27 1 and Ao and Rammer [24| each give 
approximate solutions of the influence of the envi- 
ronmental fluctuations on the Landau-Zener prob- 
ability. Using the numerically exact quasiadiabatic 
propagator path-integral (QUAPI) technique, we 
recently described the Landau-Zener probability in 
the full parameter space [2^| . For small sweep veloc- 
ities and medium to high temperatures, we have dis- 
covered non-monotonic dependencies on the sweep 
velocity, temperature, coupling strength and cut-off 
frequency which were not included in the previous 
approximate solutions. This behavior can be un- 
derstood in simple physical terms nontrivial 
competition between relaxation and Landau-Zener 
driving. 

The direct influence of environmental fluctua- 
tions on the dynamics of a driven two-state sys- 
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tem is much more pronounced when the system is 
initially prepared in the excited state since sponta- 
neous emission is even possible at zero temperature. 
Accordingly, the excitation survival probability Q 
(to end up in the excited state having started in 
the excited state), which is without environment 
strictly identical to the Landau-Zener probability, 
is strongly modified for all temperatures and bath 
coupling strengths 0, 27 . Although experimen- 



tally highly relevant, the excitation survival prob- 
ability received much less attention. In the limit 
of fast (nonadiabatic) sweeps, no influence of a 
bath is expected 0, 27 1- In the high tempera- 



ture limit under adiabatic conditions (slow sweeps), 
the two states are again driven to equal popula 



tion, Q SD = P SD = \{l - exp^A^/v)) [2J, |27j. 2. Model 



Thus, in both limits the excitation survival proba- 
bility is identical to the Landau-Zener probability, 
Q = P. However, this is quite different at inter- 
mediate and low temperatures since spontaneous 
emission drastically changes the excitation survival 
probability. Ao and Rammer [24j were the first 
to remark that in the limit of weak system-bath 
coupling, a discontinuity occurs. In the limit of 
strong coupling and low temperature, Kayanuma 
and Nakayama [27| find another simple analyti- 
cal expression, Q BC — P (l — P ) with P be- 
ing the Landau-Zener probability without environ- 
ment. Both, Kayanuma and Nakayama [27| and Ao 
and Rammer [241 ] give approximate solutions to the 
experimentally important parameter range of slow 
(adiabatic) sweeps and intermediate temperatures, 
but only weakly coupled environments. Neverthe- 
less, this parameter range is still largely unexplored. 

Beyond these direct approaches, the dissipative 
Landau-Zener problem was discussed in many more 
facets. Not being able to give a full review of 
all works, we just mention in passing that Moyer 
discussed the Landau-Zener problem for deca yin g 
states involving complex energy eigenvalues [30(. 
Extensions to three-state systems [3l[ or circuit 
QED problems [32j have been made recently. Fi- 
nally, we also mention that spin environments in 
context of Landau-Zener transitions have been dis- 
cussed by Garanin et al. [331 ] . 

In this paper, we investigate the dissipa- 
tive Landau-Zener problem in the full parameter 
range of sweep velocities, temperatures, damping 
strengths and cut-off frequency by means of the 
quasiadiabatic propagator path-integral (QUAPI) 
[3jl|3j,[3j||3j],[38( . It allows to include nonadiabatic 
as well as non-Markovian effects yielding numeri- 



cally exact results. In the next section we introduce 
the basic model. In the third section we discuss the 
time-dependent occupation probability of the two 
states during a Landau-Zener transition and in the 
following section we focus on the asymptotic popu- 
lations discussing the Landau-Zener (ground state 
survival) and the excitation survival probabilities in 
dependence of the sweep velocities, temperatures, 
damping strengths and cut-off frequency. We show 
that interesting features arise due to a competition 
between time scales associated to the Landau-Zener 
sweep and to dissipative transitions. Finally, we 
conclude with a short summary. 



A quantum mechanical two-state system which 
shows an avoided energy level crossing when driven 
is described by the Landau-Zener Hamiltonian (h = 
1) 



H L z{t) 



A 



vt 



(1) 



with the tunneling matrix element Ao and the en- 
ergy gap between the diabatic states vt, changing 
linearly in time with sweep velocity v. Here, cr XiZ 
are Pauli matrices and the diabatic states are the 
cigenstates (| \) and | 'f)) of a z . Asymptotically at 
times |i| ^> Aq/v, the diabatic states coincide with 
the momentary eigenstates of Hlz- Figure [T] plots 
the eigencnergies of Hlz (full lines) which show 
an avoided level crossing with minimal splitting Ao 
and the energies of the diabatic states (dashed lines) 
which form an exact crossing as a function of time. 

The Landau-Zener problem asks for the probabil- 
ity Po of the system to end up in the ground state 
at t = +oo, having started in the ground state at 
t = — oo (the corresponding one to end in the ex- 
cited state is given as 1 — Pn). Its exact solution 
dates back to the year 1932 [l|, 0, H 13 an d is given 

by 



Po(«,A ) = Kt (oo)|i(-oo))| 2 



1 — exp 



7rA| 
2v 



(2) 

The excitation survival probability Qq to end up 
in the excited state at t = +oo, having started in 
the excited state at t = — oo is strictly identical, 
Qo = Po for the coherent two-state problem. 

To include environmental fluctuations on 
Landau-Zener transitions, we couple Hlz diago- 
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nally to a harmonic bath [20j, |2l| , yielding 



H{t)=H LZ (t)-^J2 X ^ +b D 



k 



(3) 

with the bosonic annihilation/creation operators 
bk/b\. The bath influence is captured by the 
spectral function J(co) = 2auj cxp(— lj/u> c ), for 
which we choose here for dcfinitcncss an Ohmic 
form with the cut-off frequency w c and the cou- 
pling strength a [2(J 21 1. The Landau- Zener 



probability for the dissipative problem P = 
Tr{| f)(t |f/ooU>U {U' 1 } with the time evolu- 
tion operator Uoo = Tcxp[— i dtH{t)\ as 
well as the excitation survival probability Q = 
Tr{| I) (I \Uoo\ t)(t l^oo 1 } are now functions not 
only of Ao and v, but also of a,u c and the tem- 
perature T . In the following, we use uj c = lOAo 
unless specified otherwise. 

3. Time dependent occupation probabilities 

In this section, we explicitly consider the time- 
dependence of the population of the diabatic state 
| f) at any instant of time t, having started in | \) 
with probability one. This is given by 



P(t) = Tr{|t)(tl^U)Ul^r 1 } 



(4) 



with the time evolution operator Ut = 
Texp[— i J_ dt'H(t')]. Note that at asymptotic 
times, this quantity coincides with the standard 
Landau-Zcner probability P = P(t oo). Note 
furthermore that P(t) is the tunneling probability 
at any instant of time. Here, the dynamics of the 
quantum two-level system is described in terms of 
the time evolution of the reduced density matrix 
p(t) = TrslUtPoU^ 1 }, starting from a total initial 
density matrix po = ps <£> e~ HB ' T /Zb, where ps is 
the density operator of the quantum two-level sys- 
tem and Hb denotes the bath Hamiltonian which 
is assumed to be decoupled at to = — oo (which is 
in practice set to zero) and instantly switched on 
directly afterwards. Moreover, Zb = Tr{e _Hs / T } 
with Ub = 1- p(t) is obtained after tracing over the 
bath degrees of freedom. We calculate p(t) using 
the numerically exact quasiadiabatic propagator 
path-integral [34 , 35, 3^, 37 , 38| scheme. For details 
of the iterative tech niq ue, we refer to previous 

works in si, m, i3, m i . 

In brief, the algorithm is based on a symmet- 
ric Trotter splitting of the short-time propagator 



K.{tk+i,tk) for the full Hamiltonian into a part de- 
pending on the system Hamiltonian and a part in- 
volving the bath and the coupling parts. The short 
time propagator describes time evolution over a 
Trotter time slice At. This splitting is of course 
exact in the limit At — > but introduces a fi- 
nite Trotter error to the splitting, which has to 
be eliminated by choosing At small enough such 
that convergence is achieved. On the other hand, 
the bath degrees of freedom generates correlations 
being non-local in time. For any finite tempera- 
ture, these correlations decay exponentially fast at 
asymptotic times, thereby defining the associated 
memory time scale. QUAPI now defines an object 
called the reduced density tensor, which lives on 
this memory time window and establishes an itera- 
tion scheme in order to extract the time evolution 
of this object. Within the memory time window, 
all correlations are included exactly over the finite 
memory time r mem = KAt and can safely be ne- 
glected for times beyond r mcm . Then, the memory 
parameter K has to be increased, until convergence 
is found. Typical values, for which convergence can 
be achieved for our problem, are K < 12 and a 
reasonable choice is At ~ (0.1 — 0.2)/Ao. The 
two strategics to achieve convergence, namely de- 
creasing At and at the same time increasing the 
considered memory time r mcm = KAt, are clearly 
increasing both the needed K which results in se- 
vere demands considering that the needed computer 
power grows exponentially with K. Nevertheless 
convergent results can be obtained in a wide range 
of parameters. 

Fig.[2]shows the population P(t) vs. time for dif- 
ferent temperatures for a fixed sweep velocity v = 
0.02Aq in the weak coupling regime a = 0.0016. We 
start in the infinite past with P(t = — oo) = as 
the state | \) is fully populated, see Fig. Q] The 
Landau-Zener sweep reaches the minimal gap at 
t = 2500/Ao. Approaching this point, P(t) starts 
to increase. For larger temperatures, this increase 
is less pronounced than for lower temperatures. 
This is shown more explicitly in Fig. [31 where the 
corresponding derivative dP(t)/dt is shown. Note 
that this quantity may be viewed as tunneling rate 
for the dissipative Landau-Zcner transition. It is 
naturally more pronounced at low temperatures. 
Increasing the speed of the Landau-Zener sweep, 
the population P(t) of the | 4-) state shows some 
transient oscillatory dynamics before a stationary 
value is reached (results not shown) . We would like 
to point out that the iterative QUAPI approach 
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by construction allows the access to the full time- 
dependent Landau-Zcncr transition. In the fol- 
lowing, we focus on the stationary populations at 
asymptotic times without discussing how the sta- 
tionary state is reached in all studied parameter 
configurations. 

4. Landau- Zener and excitation survival 
probability 

In this section, we turn to the asymptotic popu- 
lations and study the Landau-Zener and the exci- 
tation survival probabilities. 

4-1. Landau-Zener probability at weak coupling 

Figure @] shows the Landau-Zener probability P 
versus sweep velocity v for different temperatures 
and for weak coupling, a = 0.0016. Clearly the 
regime with large velocities, v 3> Aq is distinguish- 
able from a regime with small velocities (adiabatic 
regime) v < Aq. At low temperatures one expects 
no sizable influence of the bath (2^, l23l [24 , I27I 
which should vanish totally at T = [25|, |26| . 
This is confirmed by our numerical results. For 
small v and low temperatures, T < Ao, we find 
P ~ 1, and thus unmodified compared to the pure 
quantum mechanical Landau-Zener result Pq (solid 
line). For increasing velocity, the Landau-Zcncr 
probability decreases rapidly and there is hardly 
any temperature effect in the considered temper- 
ature range, see Fig. |4l This observation agrees 
with results by Kayanuma and Nakayama who de- 
termined the Landau-Zener probability in the limit 
of high temperatures, Psd = |(1 — cxp(— 7rA 2 /w)) 



271 l22t l23l ] (dot-dashed line) , assuming dominance 
of phase decoherence over dissipation. For large ve- 
locities, Psd decreases as the pure Landau-Zcncr 



probability, Psd 



lz 



^ttA 2 /v, and accord- 



ingly no sizable temperature effect is expected. 

In the experimentally most relevant parameter 
range of intermediate to high temperatures, T > Aq 
and small sweep velocities, v < Aq, we find (as 
reported before [29^ a nontrivial and unexpected 
behavior of P. Besides an overall decrease of P 
with increasing temperature, we find (at fixed T) 
for decreasing velocity first a maximum of P at 
» max < Aq , then a minium at v m - m and finally again 
an increase. For decreasing temperatures, v m - ln de- 
creases, and P(v m - m ) increases. For high tempera- 
tures T > 500A , our data follow ni cely the predic- 
tions by Kayanuma and Nakayama 2jJ, [22J, [23j . 



This nonmonotonic behavior cannot be described 
in terms of pcrturbative approaches. Ao and Ram- 
mer derived temperature-dependent corrections to 
the Landau-Zener probability for low temperatures 
[li]]. They report an onset temperature T Q oc 1/v, 
above which temperature affects P. Thus, at larger 
velocities, the decrease of P due to increasing tem- 
perature starts at higher temperatures. This is in 
line with our findings of the maximum in P(v), but 
it does not account for the minimum and the sub- 
sequent increase of P for smaller v. In the limit of 
high temperatures, Psd —> 1/2 for u < Aq. Thus, 
Psd captures the decrease of P in Fig. 0] with in- 
creasing temperature, but it does not account for 
the nonmonotonic behavior for decreasing v. 

In Fig. [5j we compare our data (triangles for P 
(top) and circles for Q (bottom)) for T = 4Ao and 
a = 0.0016 with the result of Ao and Rammer (dot- 
ted lines). In fact, the latter describes qualitatively 
the maximum both for P and Q, but fails quantita- 
tively. Similarly, we were not able to match cither 
Eq. (54) of Ref. [27j] (dashed line in Fig. [5]) nor Eq. 
(40) of Ref. (Hi (dash-dash-dotted line in Fig. [5) 
with our exact data, see Fig. [5] Both describe a 
reduction of the Landau-Zener probability with in- 
creasing temperature but neither the maximum nor 
the subsequent minimum in the behavior versus v 
is predicted correctly. 

4-2. Physical picture 

The observed behavior can be understood within 
a simple physical picture realizing that the bath in- 
duces relaxation as main effect, whose time scale 
can compete with the Landau-Zener sweep veloc- 
ity. Since initially the system is in the ground-state, 
only absorption can occur, if an excitation with en- 
ergy A t = a/Aq + (vt) 2 exists in the bath spectrum 
and is thermally populated. Since A t is (slowly) 
changing with time, relaxation can only occur dur- 
ing a time window |£| < ^t r with the resonance 
time 

U = - J A 2 -A 2 , (5) 



in which the energy splitting fulfills the condition 
A t < A c = min{T, oj c } ^ as illustrated in Fig. [TJ 
In order for relaxation processes to contribute, the 
(so far unknown) relaxation time t t must be shorter 
than t r . 

For large sweep velocities, t r -C T r , relaxation is 
negligible and no influence of the bath is found as 
expected. In the opposite limit, t r 3> 7>, relaxation 
will dominate and the two levels will at any time 
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adjust their occupation to the momentary A t and 
T. Once A t > A c , relaxation stops since no spec- 
tral weight of the bath modes is available and the 
corresponding "critical" Landau-Zener probability 
can be estimated as 



splitting A r , i.e. 



P = I 

rc 2 



1 + tanh 



2T 



(6) 



For small but finite v, equilibration is retarded. The 
two levels need the finite relaxation time to adjust 
their occupation to the momentary A t and T. In 
this time, however, A t might exceed A c and then 
relaxation is not possible anymore. Thus equilib- 
rium is reached for an energy splitting in the past 
A t ' < A tc with t' < t c and t c the time when 
A tc = A c . Accordingly, P increases with decreas- 
ing v since A t changes slower and P(v — > 0) < P c , 
as observed in Fig. [U In Fig. H2 (main), we plot the 
Landau-Zener probability P(v = 0.005A§) (blue di- 
amonds) for the smallest investigated sweep veloc- 
ity and compare it with P c of Eq. (JB|) (blue full 
line) . 

Relaxation will maximally suppress the Landau- 
Zener transition when both time scales coincide, 
leading to a minimum of P at -y m ; n given by the 
condition 



ir(fmin) = T r (T , a, UJ C ) . 



(7) 



Within the resonance time window, \t\ < ^t r , 
only a single phonon absorption is likely. We 
thus can assume equilibration associated to a timc- 
averaged energy splitting A r = (2/t r ) J Q r dtA t ~ 
max{A c /2,Ao} and a resulting Landau-Zener 
probability 



P ■ = A 

1 mm n 



1 + tanh 



A r 
2T 



= P(«mm). (8) 



Subsequently, there is a maximum for a sweep ve- 
locity between v m - ln < u max < A§ . Fig. [5] (main) 
shows the QUAPI data P(v m i n ) (black circles) vs. 
T together with P m i n given in Eq. (JSj) (black dashed 
line) . 

For a fixed time and for weak coupling, 
we can estimate the decay rate out of the 
ground state using Golden Rule, T _1 (t) = 
7ra(AQ/At) cxp(— A t /a; c )n(At) with the Bose fac- 
tor n(A t ) = [exp(A t /T) - l]" 1 . For the time- 
dependent Landau-Zener problem at slow sweep ve- 
locities, we may assume that the bath sees a time- 
averaged two-level system and thus estimate the re- 
laxation rate r~ x by using the time-averaged energy 



A 2 



r r * ~ 7ra= ^~ ex P( — A r /uj c )n(A r ) . (9) 

For increasing temperature, relaxation becomes 
faster, and, accordingly, the condition for v m - mi 
Eq. 0, leading to 



2r; 



A 2 



A 2 



(10) 



is fulfilled for larger velocities. Qualitatively this 
picture describes the temperature dependence of 
v m in observed in the inset of Fig. El We plot u m i n 
(red squares) taken from the data of Fig. [Hand i> m i n 
according to Eq. ((T0|) . In detail, given the section- 
wise definition of the parameters A c and A r , the 
agreement in Fig. [5] is rather satisfactory, keeping 
in mind that there are no adjustable parameter in- 
volved. 

4-3. Excitation survival probability at weak coupling 

Figure [7] shows the excitation survival probabil- 
ity Q versus sweep velocity v for different temper- 
atures and for weak coupling, a = 0.0016. As for 
the Landau-Zener probability, the regime with large 
velocities, v 3> Aq is distinguishable from a regime 
with small velocities (adiabatic regime) v < Aq. 
For large sweep velocities, the excitation survival 
probability decreases rapidly. There is no sizable 
difference to the undamped case for all the consid- 
ered temperatures, as expected in the high temper- 
ature limit, Qsd(v ^> Aq) as well as in the low 
temperature limit for weak coupling Qq(v 3> Aq) 
and strong coupling Q sc {v ^> A§) since 

Q SC (V>A 2 ) ~ Qo(v > Ag) ~ Q S d(v > Ag) 

~ P (v » Ajj) ~ ±7r^- (11) 

At low temperatures, T < Ao, no sizable influ- 
ence of the bath is found as the data in Fig. [7] 
for T = O.OIAq almost coincide with the data 
for T = 0.5Aq (and even more so for data at 
T = 0.05A and T = 0.1A not shown in the fig- 
ure). However, for all temperatures we find the ex- 
citation survival probability to peak for sweep ve- 
locities between (0.1 — 1)Aq and a strong decrease 
for lower sweep velocities as predicted by [24|, |27| • 
For small sweep velocities relaxation of the excited 
state takes place and thus reduces the excitation 
survival probability as soon as the relaxation time 
becomes comparable or shorter than the resonance 
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time. Due to spontaneous emission relaxation is 
present at all temperatures. Furthermore we would 
expect in this picture that only at temperatures 
T > Ao the relaxation time becomes shorter due 
to induced emission by excited phonons and thus 
only for T > Ao a temperature effect should be vis- 
ible in the excitation survival probability. This is 
confirmed by the results in Fig. [7] Consistent with 
our simple physical picture introduced in the last 
section is that the sweep velocity of the sharp reduc- 
tion of the excitation survival probability coincides 
roughly with the minimum in Fig. 0] for T > Ao . 
It. as well, shifts to larger v for increasing temper- 
atures. 

Taking the sweep velocity v a = v[Q(T) = 5A0] 
(with temperature fixed) as a measure for the sweep 
velocity of the sharp reduction, we expect v a to be 
the sweep velocity for which the resonance time 
(Eq. ([5])) coincides with the relaxation time, up to 
a prefactor f a of order one which is due to the ar- 
bitrary definition of v a , i.e., 



fa.tr (V a ) = T d . 



(12) 



The inverse decay time of the excited state is now 
given as 

A 2 

t-) ~ 7ra=£ exp(-A r /a; c ) [l + n(A r )] . (13) 

Since thermal occupation is no limiting factor for 
the decay of the excited state, we have that A c = 
uj c and thus A r = cj c /2. Fig. [8] plots v a (black 
squares) versus temperature which nicely shows a 
temperature behavior ~ [l + n(A r )]. The black line 
marks the result of the best fit of Eq. (JT2J) with the 
fitting parameter f a ~ 3.04, confirming again the 
good agreement between our numerical data and 
our simple physical picture. 

In the high temperature limit, Q — > Psd (24. 27|. 
In our data, we observe that the peak in Q decreases 
with increasing temperature and Q at smallest in- 
vestigated sweep velocity increases with increasing 
temperature. Both features are in agreement with 
the expectation in the high temperature limit. For 
the smallest investigated sweep velocity, the relax- 
ation time should exceed the resonance time by 
far and thus we expect within our simple picture 
that the system fully relaxes towards the momen- 
tary equilibrium until the energy splitting exceeds 
A c = uj c . This yields the condition 



1 



tanh|'% 
2T 



(14) 



The inset of Fig. [5] compares Q c with Q(v — 

0. 005Aq) extracted from Fig. [7] and again shows 
satisfactory agreement. 

The excitation survival probability is numerically 
much harder to determine than the Landau- Zener 
probability. For example, at the sweep velocities 
v = 0.005Ao our data did not converge for tem- 
peratures below T = 4Ao and hence is not shown. 
According to our physical picture, the main effects 
of the bath occur during the resonance time win- 
dow. Outside of it, only multi-phonon processes 
can occur, which arc strongly suppressed due to the 
weak coupling. This directly affects our numerical 
scheme which is thus not sensitive to the very long 
times needed to perform a Landau-Zcncr experi- 
ment (from t = —00 to t = +00). In contrast, for 
the Landau-Zener probability all relaxation effects 
outside the resonance time window arc additionally 
suppressed by thermal occupation numbers of the 
phonons. This effect does not influence the exci- 
tation survival probability due to spontaneous de- 
cay. These multi-phonon processes outside the res- 
onance time window are, however, processes which 
need long memory in the bath and may thus be 
critical for convergence of our numerics. We should 
emphasize that our simple physical picture even ex- 
plains our stronger numerical efforts to obtain con- 
verged results for the excitation survival probabil- 
ity. 

In summary, at weak coupling our simple phys- 
ical picture explains qualitatively fully and quan- 
titatively satisfactorily both the Landau-Zener as 
well as the excitation survival probability. Both arc 
unmodified by the dissipative effects of the environ- 
ment for large velocities, v ^> A 2 ,. For slow veloc- 
ities , v < Aq, both are strongly influenced due to 
relaxation where a striking difference emerges since 
the excited state can relax by spontaneous emission 
whereas the ground state can only absorb a phonon 
when one is thermally occupied. Thus the Landau- 
Zener probability at vanishing temperature is un- 
modified by the presence of the bath whereas the 
excitation survival probability is strongly altered, 

1. e. decreases rapidly for sweep velocities smaller 
than v a determined by the competition of driving 
and relaxation. 

4-. 4- Medium to strong couplings 

Surprisingly, our simple picture still holds qual- 
itatively for stronger damping, when the Golden 
Rule is not expected to hold. Increasing a enhances 
relaxation and r r and r decrease. Thus, the 
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sweep velocity of the minimum in the Landau-Zener 
probability as well as v a in the excitation survival 
probability increase for larger coupling strengths a 
(at fixed temperature) . This is confirmed by Fig. [9] 
where P is shown for the same temperatures as in 
Fig. HI but for a larger value of a = 0.02. Fig. [TO] 
shows the corresponding excitation survival prob- 
ability for a = 0.02. The minimum of P and v a 
for Q are still observable for a = 0.02 for temper- 
atures Ao < T < 4Ao. At higher temperatures, 
only a shoulder remains for P in Fig. |H] and Q in 
Fig. [TO] does not exceed 1/2 anymore. Results for 
P and Q for even stronger coupling a = 0.2 are 
shown in Fig. [TT] For a = 0.2, the local extrema 
in P disappear, but a monotonic growth of P with 
decreasing sweep velocity is still in line with our 
simple picture. For coupling strengths a > 1/V2 
no bath influence is expected anymore 24[ for the 
Landau-Zener probability, consistent with our data. 
Therefore, we focus on a < 0.2. At the same time 
the excitation survival probability is expected to 
follow Q sc = P (l — P ) at strong coupling and 
low temperatures [27j. As seen in Fig. [TT] Q in- 
deed approaches Q sc with decreasing temperature 
for a = 0.2. 

This striking behavior of the excitation survival 
probability at strong coupling opens a road to de- 
termine the coupling strengths in experimental sys- 
tems. Under weak coupling conditions, the mini- 
mum in the Landau-Zener probability allows to ob- 
tain a. For strong coupling, the minimum vanishes 
but the excitation survival probability still shows a 
clear peak. Especially the temperature dependence 
of the peak clearly separates strong from weak cou- 
pling. For strong coupling, the peak height in- 
creases with temperature while for weak coupling, 
the peak height decreases with increasing tempera- 
ture. 

4-5. Dependence of P and Q on the system-bath 
coupling 

Next, we investigate the dependence of P and 
Q on the system-bath coupling in the regime of 
weak coupling. Fig. [TO] shows the Landau-Zener 
probability for different a at a fixed temperature 
T = 4Ao. For increasing coupling, w m i n shifts to 
larger velocities. In fact, v min depends linearly on 
a, see inset of Fig. [TO] The linear dependence is also 
predicted by our model, i.e., v mnl = 15.8oAq, in 
very good agreement with the fit v min = 17.54qAq. 
The decreasing maximum P(w max ) results from the 
shifting minimum. Another remarkable fact is that 



the Landau-Zener probability P(v m in) at the min- 
imum velocity is independent of a, as predicted 
by our physical picture, see Eq. ©. We esti- 
mate the averaged splitting A r = T/2 = 2Ao, in 
fair agreement with A r = 2.476Ao, obtained with 
P(v min ) = i[l+tanh(A r /2T)] from the data in Fig. 
1121 This agreement strongly supports our physical 
picture that relaxation dominates in the intermedi- 
ate temperature range for small sweep velocities. 

The excitation survival probability Q is shown 
in Fig. 1131 for different a at a fixed temperature 
T = 4Aq. For increasing coupling, v a shifts to 
larger velocities. As before u m in for the Landau- 
Zener probability, now v a for the excitation sur- 
vival probability depends linearly on a, see inset of 
Fig. [TO] Again, the linear dependence is predicted 
by our model, i.e., v a = 32.3oAq (taking the previ- 
ously determined factor f a = 3.04 into account), in 
very good agreement with the fit v m - m = 37.4«Aq. 
The decreasing peak height in Q results from the 
shifting v a . 



4-6. Dependence on cut-off frequency 

There is an additional time scale provided by the 
bath dynamics, which determines how fast the bath 
relaxes to its own thermal equilibrium due to the 
coupling to the system. It is given by the reor- 
ganization energy [2(| and depends on the cut-off 
frequency w c of the bath spectrum. In turn, the 
relaxation rate (O also depends on uj c and relax- 
ation is strongly suppressed when A t > u c . Fig. [14] 
shows P and Fig. [TO] shows Q for different u> c , rang- 
ing down to uo c = 0.5A. Such small values of the 
cut-off frequency describe slow bath fluctuations, a 
situation, for instance, typical for the biomolecular 
exciton dynamics in a protein-solvent environment 
[40j . With decreasing cut-off frequencies, the min- 
imum in P(v) as well as v a in Q shift to smaller 
v as qualitatively expected from Eq. © and (| 13[) 
respectively. We note that this is rather surprising 
since a small tu c also induces strong non-Markovian 
effects. At the same time, the Landau-Zener prob- 
ability P(u m in) decreases. With decreasing u> c , the 
resonance time t r and the averaged energy splitting 
A r also decrease. For cut-off frequencies T < uj c . 
we expect v m i n = 0.19 exp(— 2Aq/uj c ), in fair agree- 
ment with the fit Vmin = 0.23 exp(— 2Aq/uj c ). This 
is shown in the upper inset in Fig. 1141 The lower 
inset shows P(v m i n ) versus u c . The solid lines are 
predictions from our model and are in fair agree- 
ment with data. A similar analysis for v a was not 
possible since for small uj c the excitation survival 
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probability Q(v) did not sharply fall below | and 
thus v a could not be determined unambiguously 

5. Summary 

We have investigated the dissipative Landau- 
Zener problem by means of the numerically exact 
quasiadiabatic propagator path- integral [HI, HH [3(| 
approach for an Ohmic bath. Thereby we discussed 
the Landau-Zener probability (to end up in the 
ground state when starting in the ground state) as 
well as the excitation survival probability (to end 
up in the excited state when starting in the ex- 
cited state) . In the limits of large and small sweep 
velocities and low temperatures, our results coin- 
cide with analytical predictions [24, 27, 25, 32j ■ In 
the intermediate regime, when the sweep velocities 
are comparable to the minimal Landau-Zener gap 
and intermediate temperatures, we have identified 
novel non-monotonic dependencies of the Landau- 
Zener probabilities on the sweep velocity, temper- 
ature, system-bath coupling strength and cut-off 
frequency. This parameter range is clearly not ac- 
cessible by perturbative means. The observed be- 
havior can be understood in rather simple physi- 
cal terms MS i\ nontrivial competition between re- 
laxation and Landau-Zener driving. The main dif- 
ference between the Landau-Zener and the excita- 
tion survival probability results from the simple fact 
that the excited state can always decay via spon- 
taneous emission while the ground state needs a 
thermally excited phonon in order to become ex- 
cited. Thus even at vanishing temperature the ex- 
citation survival probability decreases strongly for 
small enough driving speed whereas the Landau- 
Zener probability only for high temperatures. As 
nowadays advanced experimental set-ups allow for 
a rather comprehensive control of the parameters, 
this novel feature should be accessible by available 
experimental techniques. 

We thank V. Peano and S. Ludwig for discussions 
and acknowledge support by the Excellence Initia- 
tive of the German Federal and State Governments. 
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in time through an avoided crossing. The solid lines mark 
the exact energies whereas the dashed lines give the crossing 

energies ±vt of the diabatic states. In addition, the reso- Fi § ure 3: Thc timc derivative dP/dt corresponding to Fig. 

nance time window of width t r around t = is shown, in El 
which At < A c , see Sec. 14.21 for more details. 




Figure 2: The occupation probability P(t) of the diabatic 
state | f) versus time and temperature is shown for a slow 
sweep velocity v = 0.02Aq, a weak system bath coupling 
a = 0.0016 and a large cut-off frequency ui c = 10 Ao. The 
time is restricted around the crossover which was at t CIOBB = 
2500/ A . 



Figure 4: The Landau-Zencr probability P for various tem- 
peratures T is shown for a weak system-bath coupling a = 
0.0016. The dotted lines are guides to the eye. The solid line 
marks the coherent Landau-Zcner probability Po, while the 
dot-dashed line indicates the high-temperature limit Psd> 
see text. 
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Figure 5: Comparison of the Landau-Zcner (P, top panel, 
triangles) and the excitation survival (Q, bottom panel, cir- 
cles) probability with perturbative approaches for T = 4Ao 
for a weak system-bath coupling a = 0.0016. The thick 
solid line in both panels marks the standard Landau- Zener 
probability Po, while the dot-dashed line shows Pso- The 
dotted lines show the behavior according to Ao and Rammer 
241] , the dashed lines according to Kayanuma and Nakayama 
2711 . Moreover, in the top panel, the dash-dash-dotted line 
marks the result of Pokrovsky and Sun 28] (which only give 
a formula for the Landau- Zener probability P). In the lower 
panel, the dot-dot-dashed line indicates Pq(1 — Po). 




Figure 7: The excitation survival probability Q for various 
temperatures T is shown for a weak system-bath coupling 
a = 0.0016. The dotted lines are guides to the eye. The solid 
line marks the coherent excitation survival probability Qo = 
Po which is identical to the Landau-Zcner probability. The 
dot-dashed line indicates the high-temperature limit Psd 
and the dot-dot-dashed line indicates the strong coupling 
limit Q BC = P (l - P ). 
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Figure 6: Main: Landau-Zcner probability P(v = 0.005Aq) 
at the smallest sweep velocity (blue diamonds), P c of Eq. ((5J 
(blue solid line), Landau-Zener probability P(f m i n ) at the 
sweep velocity of the minimum (black circles), and P m in °f 
Eq. ((Sj (back solid line) versus temperature. Inset: The 
minimum sweep velocity (red squares) extracted from the 
data of Fig. [4] and n m i n according to Eq. 1101 (red dotted 
line). The coupling strength is a = 0.0016. 
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Figure 8: Main: Sweep velocity v a = v[Q(T) = ^Ao] (with 
temperature fixed) for which the resonance time (Eq. JSJ) 
coincides with the relaxation time (up to a prefactor f a ) vs. 
temperature. Inset: Critical excitation survival probability 
Q c , Eq. flit (solid line), and Q(v = 0.005Ag) extracted from 
Fig. [7] versus temperature for a = 0.0016. 
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Figure 9: The Landau-Zcncr probability P for various tem- 
peratures T is shown for an intermediate value of the system- 
bath coupling a = 0.02. The dotted lines are guides to 
the eye. The solid line marks the coherent Landau-Zener 
probability Po, while the dot-dashed line indicates the high- 
temperature limit Psd- 



Figure 11: The Landau-Zener (P, all data above the red 
dash-dotted line) and the excitation survival (Q, all data be- 
low the red dash-dotted line) probability for various temper- 
atures T are shown for strong system-bath coupling a = 0.2. 
The dotted lines are guides to the eye. The solid line marks 
the coherent excitation probability which is identical to the 
coherent Landau-Zener probability Pq. 
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Figure 10: The excitation survival probability Q correspond- 
ing to the case shown in Fig. [9] Here, the dot-dot-dashed 
line marks the strong coupling limit Q sc = Po(l — Po). 
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Figure 12: Landau-Zener probability P for different a for 
T = 4Aq. Inset: Linear dependence of u m i n on a, see text. 
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Figure 13: Excitation survival probability Q for different a 
for T = 4Aq. Inset: Linear dependence of v a on a, see text. 
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Figure 15: Excitation survival probability Q corresponding 
to Fig. [H] 



Sweep velocity v / (A ) 



Figure 14: Landau-Zcncr probability P for various cut-off 
frequencies u) c for T = 4Ao and a = 0.00985. Insets: 
^(^min) (bottom) and u m in (top) vs u) c \ solid lines: model 
predictions, see text. 
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